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A theory of the proximity effects of the exchange splitting in a ferromagnetic metal on a two 
dimensional electron gas (2DEG) in a semiconductor is presented. The resulting spin-dependent 
energy and lifetime in the 2DEG create a marked spin-splitting in the driven in-plane current. The 
theory of the planar transport allows for current leakage into the ferromagnetic layer through the 
interface, which leads to a competition between drift and diffusion. The spin-dependent in-plane 
conductivity of the 2DEG may be exploited to provide a new paradigm for spintronics devices 
based on planar devices in a field-effect transistor configuration. An illustrative example is provided 
through the transport theory of a proposed spin-valve which consists of a field-effect transistor 
configuration with two ferromagnetic gates. Results are provided for two experimentally accessible 
systems: the silicon inversion layer and the naturally-formed InAs accumulation layer. 
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I. INTRODUCTION 



The emerging field of spintronics aims to implement 
semiconductor devices which utilize both the carrier 
charge and spin degrees of freedomi Research in the 
field has been inspired by the early device proposal of 
Datta and Das,^ which consists of spin injection through 
a ferromagnetic-semiconductor interface and spin ma- 
nipulation using the Rashba spin-orbit effect There 
has been recent progress in achieving maximum spin 
injectio n^i^i^i^ and in characterization of the Rashba 
effec t^i?'^*^'^^ . However, a spin device based on injection 
has a high tunnel resistance, analogous to a Schottky 
diode. We have suggested an alternate approach of spin 
creation and manipulation which makes use of the prox- 
imity effects of a ferromagnet on a semiconductor li^ The 
eventual device would avoid the high tunnelling resis- 
tance by keeping the main driven current path entirely 
in the semiconductor through normal ohmic nonmagnetic 
leads, resembling the field-effect transistor design. 

In this paper we present a comprehensive theory of the 
consequences of the ferromagnetic proximity on the equi- 
librium and transport properties of the electrons in the 
semiconductor. First, we fully examine the coupling of 
a semiconductor two dimensional electron gas (2DEG) 
with a ferromagnetic layer through a very thin poten- 
tial barrier. Explicit calculations are provided for two 
realizable systems, the silicon inversion layer and the 
naturally-formed InAs surface layer, both with ferromag- 
netic gates. The coupling is conveniently treated by a 
Green's function method which can account for realis- 
tic confinement for the semiconductor 2DEG as well as 
the effects of a short electron mean free path in the fer- 
romagnet. The complex self-energy of the 2DEG due to 
the interaction with the ferromagnet contains both static 
and dynamic effects (a Zeeman-like splitting and spin- 
dependent scattering times, respectively). Both proper- 
ties alter the in-plane conductivity and can be exploited 
through the spin-dependent transport under the ferro- 
magnetic gate. As an illustration of the possible con- 



sequences of the ferromagnetic proximity on the semi- 
conductor, we show in detail the behavior of the density 
and current in a previously proposed^^ planar spin-valve 
in MOSFET-style configuration with two adjacent ferro- 
magnetic gates with reversible magnetizations. We find 
that the proposed device has a significant magnetoresis- 
tance effect for reasonable system parameters. 

The coupling of the 2DEG to the ferromagnetic gate 
exponentially decreases with their separation by a poten- 
tial barrier. Our study includes oxide barriers down to 
the smallest state-of-the-art thickness and direct contact 
between the 2DEG and the ferromagnet to ensure sizable 
spin effects. A consequence of this requirement for ultra- 
thin oxides is that current leakage into the ferromagnetic 
gate(s) becomes an integral part of the in-plane electron 
transport. This naturally results in the need to consider 
both the current driven by the in-plane electric field and 
the diffusion driven by the density gradient, analogous to 
the bulk case >"? ' '^" i Although such leakage effects would 
decrease the efficiency of normal field-effect transistors, 
we find that they actually enhance the spin-dependence 
of the electron transport from the source to the drain. 

We argue that the relevance of proximity effects be- 
tween a semiconductor and a ferromagnet is supported 
by recent experiments. A series of optical Faraday rota- 
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have shown that unpolarized non- 



equilibrium electrons in a semiconductor can sponta- 
neously acquire a net spin polarization in the presence 
of a ferromagnetic interface. We have interpretedii the 
observed effect as arising from the spin-dependent reflec- 
tion of electrons off the ferromagnetic interface. Assum- 
ing negligible spin scattering at the interface, the spon- 
taneous spin polarization produced in the semiconductor 
is indicative of the strength of the coupling across the 
semiconductor-ferromagnet junction. In addition, replac- 
ing the Schottky barrier-^'-^ with a much thinner oxide 
will only increase the interaction. 

Proximity effects between dissimilar materials have 
been used to describe how the ordered state of one 
medium induces a similar order in the other medium 
which is otherwise normal. The induced order parameter 



decays away from the interface. The common examples 
are the proximity effects between a superconductor and 
a normal metal (or a semiconductor 2DEG) and between 
a ferromagnet and a paramagnetic metal. The induced 
ferromagnetic order in a non-magnetic metal is weak and 
requires the metal to be superparamagnetic^ such as Pd 
or Pt. In this paper we concentrate on the proximity ef- 
fect between a ferromagnet and a semiconductor 2DEG. 
The 2DEG is examined rather than the bulk semicon- 
ductor because it has two advantages: (1) the 2DEG 
is confined near the interface where it is more suscep- 
tible to the influence of the exchange-split electrons in 
the ferromagnet, and (2) while the induced polarization 
in the 2DEG may be too small for a magnetization mea- 
surement, the influence on the spin-polarized transport is 
the ultimate goal. Moreover, there exists a large amount 
of knowledge and technology dedicated to the manipu- 
lation of semiconductor 2DEG's, further increasing the 
potential for spintronics devices. 

In Sec. mwe explain the model of the semiconductor 
2DEG-ferromagnet used in the calculations, followed by 
the derivation of an effective tight-binding Hamiltonian 
(Sec. Ill All which is used in a Green's function formal- 
ism (Sec. Ill B|) to approximate the coupling between the 
semiconductor 2DEG and the ferromagnet. In Sec. IIIII 
we calculate the self-energy of the 2DEG coupled to the 
ferromagnet for two experimentally realizable systems, 
the silicon inversion layer (Sec. IIII Ap and the naturally- 
formed InAs surface layer (Sec. IIIIIj|l . In Sec. IIVI we 
derive the equations and show representative results for 
the 2DEG density and in-plane current for the silicon in- 



version layer (Sec. IIV All and for the InAs surface layer 
(Sec. HV^ . In Sec. |V|we calculate the 2DEG density 
and current for a spin- valve with two ferromagnetic gates 
for the sihcon (Sec. IV A|l and the InAs (Sec. IVB|I sys- 
tems. We summarize our findings and draw conclusions 
in Sec. EH 



II. MODEL OF THE ELECTRONIC 
STRUCTURE 

The two particular systems we will consider have band 
diagrams as shown in Fig. ^ The semiconductor 2DEG 
is on the right, which is induced by a gate bias for the 
silicon inversion layer (Fig.^a)) and is naturally formed 
for the InAs surface layer (Fig.^b)). The left side is the 
ferromagnet, which is modelled as parabolic bands split 
by the exchange energy A.^^ A very thin oxide separates 
the 2DEG from the ferromagnet. This oxide is necessary 
in the silicon system because a gate bias is necessary for 
the creation of the inversion layer. Because the surface 
layer in the InAs system forms naturally, no gate bias 
is necessary and hence intimate contact (no oxide bar- 
rier) is possible between the InAs surface layer and the 
ferromagnet. 

All calculations are done within the effective mass 
approximation. This, along with the exchange-split 
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FIG. 1: Band diagrams for (a) the silicon inversion layer and 
(b) the InAs surface layer, both separated from a ferromag- 
net by a thin oxide barrier. The ferromagnet, with exchange 
splitting A, is on the left and the semiconductor 2DEG is on 
the right. In (a), the Fermi level in the ferromagnet is lower 
than the Fermi level in the silicon inversion layer. In (b), the 
Fermi level in the ferromagnet is equal to the Fermi level in 
the InAs surface layer. 



parabolic band approximation, cannot possibly account 
for all band effects in the ferromagnet. We wish to show 
in a transparent way how the spin-dependence in the fer- 
romagnet can influence the 2DEG; more realistic calcula- 
tions will be necessary for comparison with experiment, 
but the essential effects we deduce here should remain 
valid. In addition, the two spin channels are considered 
to be completely decoupled and are labelled by + and 
— . Experimentally, it appears that the spin lifetime in 
semiconductor heterostructures can persist for long times 
and over long distances, including through heterostruc- 
ture interfacesj^SiSiiSSiS^ We neglect any mixing due to 
the spin-orbit effect or other spin-flip processes. The ne- 
glect of the Rashba spin-orbit interaction^ is valid in Si 
and holds less well in the InAs surface layer, but suffi- 
ciently well without gate voltagesS 

In a recent articlei^ we reported the calculation of the 
coupling between a silicon inversion layer and a ferromag- 
netic gate in the effective mass approximation using a tri- 
angular potential in the semiconductor region. The sim- 
plicity of the potential allowed us to use the wavefunction 
matching conditions to derive an equation specifying the 
complex energy of the eigenstate of the coupled system. 
Strong scattering in the ferromagnet was incorporated by 
putting a phenomenological damping into the wavefunc- 
tion in the ferromagnet region. The triangular potential 
approximation, although valid in the silicon system for 
strong inversion, is not appropriate for the weakly con- 
fined 2DEG at the InAs surface. For this reason, the InAs 
potential is calculated self-consisently with the surface 
layer density distribution using the coupled Schrodinger 
and Poisson equations. A Green's function method is 



devoloped to calculate the coupling between the semi- 
conductor 2DEG and the ferromagnet which works for 
both the triangle potential in the silicon system and the 
self-consistent potential in the InAs system. In addition, 
the Green's function approach is ideal because it allows 
us to include scattering in a more natural and rigorous 
manner than in the wave function approach. In brief, 
the coupled semiconductor 2DEG-ferromagnet system is 
solved in three steps. (1) Separate the semiconductor 
and ferromagnet regions and solve them separately. (2) 
Approximate the coupling between the two subsystems 
by transforming the original effective mass Hamiltonian 
into a "tight-binding" (i.e., tunneling) Hamiltonian. (3) 
Calculate the self-energy of 2DEG electrons due to the 
interaction with the ferromagnet. 



A. Derivation of the Effective Hamiltonian 

The effective mass Hamiltonian for the ferromagnet- 
oxide-semiconductor junction (see Fig. 
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FIG. 2: A cartoon of the splitting up of the potential. The 
original band diagram is on the left, and this is split into a 
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Fermionic anti-commutation relations between the oper- 
ators in all parts of the system, we orthogonalize the fer- 
romagnetic eigenstates to the discrete 2DEG eigenstates. 
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+Ubiz)e{zy, - z)e{z) + UscizMz - z^) , (i) 

can be written as H ~ K + V'^'^ + Vj™, where K is the 
kinetic energy and the semiconductor and ferromagnet 
potentials are defined respectively as 

v^'^iz) = c/,c(z)e(z-Zb) + o-e(zb-z) , 
vj"'{z) = [/b(z)e(zb-z)e(z) + o-e(z-Zb) 
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cr -M e(-z) 
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The zero of energy is put at the right side of the barrier 
for convenience. Basically, V"J'"(z) consists of the spin- 
dependent ferromagnet potential (c/fm + Y'^-M) for 

z < 0, the oxide barrier potential Uh{z) for < z < Zb, 
and zero for z > Zb. From now on we choose the spin 
quantization axis to be parallel to the magnetization M 
in the ferromagnet, so that V^™ The semicon- 

ductor potential V^'^{z) consists of the 2DEG confine- 
ment potential [4c (z) for z > Zb and is zero everywhere 
else. 

This separation of the potential (illustrated in Fig. ^ 
is convenient because it allows us to solve the two subsys- 
tems with potentials y^™(z) and V^'^{z) separately, and 
then calculate the coupling between them. The subsys- 
tems are solved, yielding the semiconductor eigenstates 
{Xni^)} with energy and the ferromagnet eigenstates 
|x|™j_(z)| with energy e[™_j_. The finite overlap between 

the wavefunctions for the two subsystems implies that 
they are not mutually orthogonal. To ensure proper 



This allows us to write a general state of the system as 



k 



(4) 



where a„ ± and Ck^± are the Fermion operators for the 
2DEG state X^i^) ^-nd the orthogonalized ferromagnet 
state (f>^^j. (z) respectively. The only Fermion operators 
that do not anticommute properly are 



'Jk,k' ■ 



(5) 



As an approximation we keep only terms that have a 
single exponentially decaying term in an overlap integral. 
Since {Xn\x^k^±) iXk^ ±\Xn) contains the product of two 
overlap integrals with an exponentially decaying term, 
we neglect it, leaving an approximate set of Fermionic 
operators. 

Now we must calculate all the matrix elements involved 
in {tp±\K + VjI^ + V^''\tp±), keeping the approximation as 
mentioned above. The matrix elements between different 
semiconductor states are 



{xl',\K + Vl^^ + VnXn) 



fScsj I yfm 



where 



vx^^= I dzixl^rV^X 



fm. 



(6) 



(7) 



is neglected because the integral contains two exponen- 
tially decaying functions. The matrix elements between 
a semiconductor state and a ferromagnet state are 



^^\K + Vt^ + V^'^\xi:±) 



bt%\xl 



(8) 
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where {(j)^™j_\Xn) = because of the orthogonahty be- 
tween the semiconductor and ferromagnet states, and 

C±= / dzicl^'^^Yvt'xl' (9) 

is kept because there is only a single exponentially decay- 
ing wavefunction in the nonzero integration region. Fi- 
nally, the matrix elements between ferromagnetic states 
are 

« e'^Ju,k' (10) 

where e^^-j-is the energy of the state Xfc"± before the or- 
thogonalization to the 2DEG states and 

n 

+ (C(X^„1X1!±) + K.'",±) {X'^,±\x'n)] (11) 

is neglected because all terms contain more than a single 
exponentially decaying function. 

We are left with an effective Hamiltonian with a simple 
tight-binding form: 

n,(7 k,a 

+ J2 {<.Vnl^Ck,a + clXl.^n.^) (12) 
n,/c,cr 

This Hamiltonian represents the coupling between the 
quantum confined electrons in the semiconductor and the 
broad spin-dependent continuum in the ferromagnet. We 
solve the Hamiltonian below using a Green's function ap- 
proach. 

B. Green's functions for the Coupled System 

We are interested in the properties of the 2DEG elec- 
trons due to the coupling with the ferromagnet. The 
retarded Green's functions relevant to our calculations 
are: 

Gn,±{t) = -*e(t)({a„,±(t),al_±(0)}) , (13) 

Gfe,±(t) = -»e(t)({cfe,±(t),c^(0)}) , (14) 

Gk,n,±{t) - -*e(i)({cfe,±(t),a,\,±(0)}) , (15) 

where Q(t) ensures us that t > 0. The unperturbed 
Green's functions when the coupling V^™ _j_ is zero are 

Gi%(i?) = {E-e:: + ^0+y' , (16) 
G^^^iE) = {E-e'^X+n',:%y' ■ (17) 



Note that the ferromagnet is assumed to be a dirty con- 
ductor, so that an imaginary part has been added to the 
ferromagnet energy to account for strong spin-dependent 
scattering in the ferromagnet ef™j_ — > e|"_j_ — «7^°± . 

Calculating the Green's functions with the interaction 
from the Hamiltonian and Fourier transforming, we have 
the following equation for the full 2DEG Green's func- 
tion: 

G„,± (E) = {E~ el' + tO+ - I]„,± {E)) , (18) 
where the self-energy of the 2DEG electrons is 

^n,±{E) = E Vnl±Gt±iE)Vkl± ■ (19) 

k 

The effect of the coupling on the 2DEG eigenstates is 
found by from the complex energy E which satisfies 

^-e^,^ + ^0+-S„,±(^) = . (20) 

The real part of the 2DEG self-energy, A„^-|-(£') — 

Ke (^T,n.±{E)j , is the spin-dependent level shift of the 

2DEG subbands due to the coupling with the ferro- 
magnet. This gives the 2DEG a static spin-splitting, 
which is |A+ — A_| if only one 2DEG subband is occu- 
pied. The two spin channels will have different densities, 
^ N_ . The imaginary part of the 2DEG self-energy, 

h/2Tn^±{E) = — Im (T,n,±{E)j is related to the lifetime 

of the 2DEG electrons to scatter in a spin-dependent way 
off the ferromagnet. This new spin-dependent scattering 
channel will result in different conductivities for the two 
spin channels, and will be addressed in Sec. IIVI 

III. RESULTS FOR RELEVANT SYSTEMS 

We have shown how to calculate the spin-dependent 
effects for the coupled 2DEG-ferromagnet system. In 
this section we apply this method to two experimentally 
realizable systems: the silicon inversion layer and the 
naturally-formed InAs surface layer. Silicon MOSFET- 
type devices are ubiquitous in present technology and any 
spintronics device in this system would have a strong in- 
dustrial base. The inversion layer is created by a gate 
bias, which increases the coupling with the ferromagnet 
by pressing the electrons close to the interface with the 
(ferromagnetic) gate. The need for ultra-thin oxide bar- 
riers (to increase the coupling) implies that 2DEG elec- 
trons can irreversibly tunnel into the ferromagnet. We 
shall account for the effects of the current leakage. Since 
the removal of the oxide barrier (provided that it is not 
replaced by a Schottky barrier) brings simplicity, we also 
examine the naturally occurring InAs accumulation layer 
which forms an ohmic contact with the metal. The 2DEG 
forms at the interface without the need for a gate bias, so 
the leakage of electrons is not a problem in this system. 
However, the confinement of electrons in such a surface 
layer is quite weak. 
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FIG. 3: The ferromagnet subsystem, approximated as a large 
but finite box. It consists of the exchange split ferromagnet 
potential, the barrier potential, and is zero for z > z^. The 
wavefunction must vanish at the left boundary z = — ^fm• 
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We start with calculation of the ferromagnet subsys- 
tem wavef unctions and energy eigenvalues; we will show 
below how to treat the semiconductor subsystem. The 
ferromagnet subsystem is made of the ferromagnet po- 
tential for z < 0, the barrier potential for Q < z < z^, 
and is zero for z > Zh- AH states are included in the cal- 
culation that decay exponentially for z > Zb ( all states 
with energy less that zero). To fix the normalization of 
the ferromagnet states and the ferromagnet density of 
states, the ferromagnet continuum is approximated by a 
large but finite box (see Fig. |3Jl. Convergence is checked 
for both the size of the box and the number of points in 
the box. 

The ferromagnet has the same parameters for all calcu- 
lations below. The exchange-split parabolic bands have 
an effective mass equal to the free electron mass mo. The 
majority wavevector is k!^]^ = 1.1 A and the minority 

wavevector is k^I^ — 0.42 A ^ ^ corresponding to a ma- 
jority Fermi energy of C/fm = 4.6 eV and an exchange 
energy of A = 3.9 eV. The ferromagnet is assumed to be 
dirty, so that the scattering is significant. We account for 
this by putting by hand an imaginary part into the ferro- 
magnet energy eigenvalues after calculating the particle- 
in-a-box eigenstates that is the same for all wavevectors, 
7*™ = 1.1 eV and 7!"° = 0.8 eV,^ We are left with 
the wavefunctions and complex energy eigenvalues for all 
states in the ferromagnet with the real part of their en- 
ergy less than zero. Below we describe separately how 
we calculate the semiconductor subsystem for the two 
systems we consider in this paper. 



A. Silicon Inversion Layer 

The 2DEG in this system is strongly inverted at the in- 
terface and the triangular approximation to the semicon- 
ductor potential is adequatei^ The dielectric constant in 
silicon is taken as esi = 11.7, the longitudinal effective 
mass responsible for the confinement is rn^-^y = .91 mo, 
and the transverse effective mass relevant to the in-plane 
motion of the 2DEG is mgj ^ = .19 mp. The oxide barrier 



FIG. 4: A plot of the spin splitting | A+ — A_ | (dotted line run- 
ning from upper left to lower right) and the spin-dependent 
scattering times t+ (full line running from lower left to upper 
right) and r_ (dashed line running from lower left to upper 
right) as a function of the thickness of the Si02 barrier for a 
silicon inversion layer. The parameters used are explained in 
the text. 



between the silicon and the ferromagnet is assumed to 
be Si02 with barrier height (from the ferromagnet Fermi 
level) ?7si02 ~ 3.2 eV, effective mass "igjQ^ = .3 mo, 
and dielectric constant £3102 = 3.9i^ We assume that 
the electric field in the barrier, which is set by the bias 
between the silicon substrate and the ferromagnetic gate, 
is near its breakdown value of £'si02 — 12 MV/cm. The 
electric field in the silicon responsible for the 2DEG con- 
finement is then E'gi = (esi02/£si) ■ -E'Si02- The density in 
the 2DEG is assumed to be 10^^ cm~^, corresponding to 
a Fermi energy of w 6.3 meV. 

The numerical results for the spin-splitting |A+ — A_| 
(dotted line running from upper left to lower right) and 
the scattering times for the two spin channels T4. (full line 
running from lower left to upper right) and t_ (dashed 
line running from lower left to upper right) are shown 
in Fig. ^for -Esi02 = 10 MV/cm. The horizontal axis 
is the thickness of the oxide barrier; sizable coupling be- 
tween the 2DEG and the ferromagnet occurs only for 
ultrathin oxides at the limit of current device fabrication 
techniques. The spin-splitting is very small even for the 
thinnest oxide considered; at 6 A barrier width, the spin- 
splitting of « 0.1 meV versus the 2DEG Fermi energy of 
w 6.3 meV yields a static spin polarization of the 2DEG 
« 1%. 

On the other hand, the scattering time for the two spin 
channels is a promising effect of the coupling that could 
be used in a device. For 6 A barrier width, the two spin 
channels have scattering times of approximately 3 ps and 
6 ps; the spin-dependent coupling to the ferromagnet has 
opened a new spin-dependent scattering channel for the 
2DEG electrons. To utilize this effect in a device, the 
new scattering times must be comparable to the intrinsic 
scattering time for the silicon inversion layer (including 
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phonons, impurities, defects, etc., but not including the 
effects of the nearby ferromagnetic layer), which at low 
temperatures is near 1 psf^ As is evident in the plot, 
the scattering times for the two spin channels approach 
this intrinsic value for very thin oxide thicknesses, while 
their relative ratio remains the same. Because the two 
times are near the intrinsic time, but are still quite differ- 
ent (more than a factor of 2), we would expect that any 
physical quantity that depends on the scattering time 
would see this spin effect. This will be addressed later 
when we discuss 2DEG transport under the gate. 



B. InAs Surface Layer 

The 2DEG at the InAs surface is naturally formed. 
As opposed to the strongly inverted silicon inversion 
layer discussed above, the confining potential in the InAs 
system is quite weak and a realistic calculation of the 
coupling requires that the potential be calculated self- 
consistently with the density distribution in the 2DEG. 
We use the coupled Schrodinger and Poisson equations 
in the semiconductor region z > Zh, 



> 

& -1.3- 



n'^ d ( I dxn 



2 dz 



dz 
92 y: 



y-(z)x^=(z) = e-X„^^(zX21) 



dz^ 



-NoJ2\x;.%z)\^ (22) 



where the total density in the 2DEG is kept fixed, A'q = 
(2K„As/'^^^)E„(eF - '^n)- The band bending in the 
barrier is negligible. Unlike the silicon inversion layer, 
more than one subband in the InAs surface layer is oc- 
cupied by electrons. Solving these equations simultane- 
ously results in the 2DEG wavefunctions X^i^)^ energies 
e^'^, and the semiconductor potential V^^{z). The InAs 
parameters used in the calculation are an effective mass 
of 

"^inAs = -023 mo, dielectric constant einAs = 14.6, 
and 2DEG density A^o = 10^^ cm'^. The barrier is 
taken as AI2O3 with a barrier height from the ferro- 
magnet Fermi level of C/AI2O3 — 1-2 eV, effective mass-- 



TO 
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.75 Too, and dielectric constant eAbOa = 3.9. 



The wavefunctions and energies associated with the 
calculation of the 2DEG with a 10 A barrier are shown 
in Fig. [S] As is evident in Fig. EJa), the confining 
potential for the 2DEG is quite shallow and weak; the 
depth of the potential is only a few hundred meV and the 
potential does not flatten out until about 300 A from the 
interface. Two subbands are occupied up to the Fermi 
level. The wavefunctions of the two subbands are shown 
in Fig. Efb), and the density distribution in the 2DEG is 
shown in Fig. [Sfc). 

The numerical results for the coupling of the first sub- 
band of the InAs surface layer with a ferromagnetic gate 
are shown in Fig. as a function of the AI2O3 barrier 
thickness. The coupling for the second subband is an 
order or magnitude weaker, and because only « 20% of 
the carriers are in this subband, we neglect its presence 
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FIG. 5: (a) The calculated potential (eV) for the InAs sur- 
face layer with the parameters given in the text. The zero of 
energy is at the top right side of the barrier, (b) The wave- 
functions for the two subbands (arb. units), (c) The density 
distribution in the 2DEG (arb. units). 




FIG. 6: A plot of the spin splitting | A+ — A_ | (dotted line run- 
ning from upper left to lower right) and the spin-dependent 
scattering times t+ (full line running from lower left to upper 
right) and r_ (dashed line running from lower left to upper 
right) as a function of the thickness of the AI2O3 barrier for 
the lowest subband in the InAs surface layer. The parameters 
used are explained in the text. 



from now on. The spin splitting is negligible except for 
intimate contact, in which case the value of « 10 meV 
means a polarization of the gas of ~ 15%. This value 
is quite high and is achievable in this system because 
the barrier is unnecessary for the creation of the 2DEG. 
The scattering times for the two spin channels approach 
the sub-picosecond range for the thinnest barriers and 
intimate contact. The intrinsic scattering time of the 
surface InAs 2DEG is in the hundreds of fs range at low 
temperaturesi^^ The coupling has been calculated down 
to an oxide thickness of 1 A; the tight-binding model 
breaks down for intimate contact. 
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IV. SPIN-DEPENDENT TRANSPORT UNDER 
THE FERROMAGNET 

We have calculated the spin-dependent self-energy for 
2DEG states coupled to a ferromagnetic gate, which in 
general results in (1) a spin-dependent scattering time as- 
sociated with the interaction of the 2DEG electrons with 
the ferromagnet t± and (2) a spin splitting ( |A+ — A_| 
in the single subband limit) resulting in unequal electron 
densities in the two spin channels iV-|- . In this section we 
examine the influence of these two results of the coupling 
on the in-plane transport of the 2DEG, and find that (1) 
the in-plane conductivity becomes spin-dependent and 
(2) 2DEG electrons can irreversibly leak into the ferro- 
magnetic gate in a spin-dependent manner (if the ferro- 
magnet is biased with respect to the semiconductor). 

We derive the transport equations in the 2DEG, ac- 
counting for current leakage into the gate and, therefore, 
the density variation along the semiconductor channel. 
In the following we assume that the current in the 2DEG 
flows in the ^-direction. The boundary condition on the 
channel density will be that at the source and drain con- 
tacts the 2DEG takes its equilibrium density related to 
the confinement and the Fermi energy of the 2DEG. The 
growth axis is in the i-direction, with the ferromagnet 
interface at z = 0. The electron confinement in the z- 
direction is assumed to be constant along the channel. 
The system is homogeneous in the j/-direction, so there 
is no y dependence in any of the equations. The two 
spin channels are considered to be completely decoupled 
throughout the device (valid for channel lengths smaller 
than the spin relaxation length). 

The continuity equation for 2DEG electrons is 

d - ^ 

—n±{x,z,t) + V-j±{x,z,t) = 0, (23) 

where n±{x,z,t) is the spin-dependent particle density 
( cm~'^) and j± {x,z,t) is the spin-dependent particle cur- 
rent (cm~^ ■ s^^). In the steady-state, 



-T^JxMx, z) = -Tr]z.±{x, z) 
ax oz 



(24) 



Integrating out the z-dependence through the semicon- 
ductor up to the interface at z = we obtain 



d_ 

dx 



Jx,±{x) = -3z,±{x,Z ^Q) + j:,^±{x,Z = -QO) 



-iz,±{x,z = 0) 



(25) 



where Jx,±{x) = J^^dzjx,±{x, z,) is the integrated 
2DEG current flowing in the ^-direction. The term 
jz,±{x,z = 0) represents the leakage of 2DEG electrons 
irreversibly into the ferromagnetic gate and will be dis- 
cussed separately for the silicon and InAs systems below. 
The term jz,±{x,z = — oo) represents the electrons that 
are injected from the semiconductor substrate into the 
2DEG. This process is unfavorable because electrons are 



injected into the 2DEG much more efficiently from the 
source and drain contacts, so that jz±{x, z = — oo) = 0. 
The Drude conductivity is ordinarily 



(26) 



where A'o is the 2DEG density (cm~^), — e is the electron 
charge, tq is the intrinsic lifetime, and m* is the effec- 
tive mass of the electrons. The coupling of the 2DEG 
with the ferromagnet introduces spin-dependence to the 
2DEG density Nq and lifetime tq, and, hence, to the con- 
ductivity. 



iV±e^f± 
m* 



(27) 



where N-i-{x) — j'^^dzn±{x, z). The spin-dependent 
lifetime f± includes both the intrinsic 2DEG lifetime tq 
and the spin-dependent scattering time associated with 
the ferromagnet t±, 



1 

To 



1 



(28) 



The spin dependence in the conductivity will affect the 
transport of the 2DEG only if the spin-dependent den- 
sities N± are very different or the spin-dependent scat- 
tering times T± are very different and are comparable to 
the intrinsic 2DEG lifetime tq. 

Using this conductivity, the in-plane charge current is 

-eJ,,±(a;) = a±{x)E,, + eD±{x)-^N±{x) . (29) 

The diffusion constant D±{x) can be related to the con- 
ductivity using the Einstein relation, and at finite tem- 
perature in two dimensions is 



D±{x) = 



[{x)f± 
m* 



1 



,-4(x)/fcBT 



(30) 



where f^{x) — [4:nN±{x)] is the chemical potential 
at T = 0. The in-plane field. Ex, results form a source- 
drain bias and is approximately constant throughtout the 
2DEG (we do not consider the feedback of the 2DEG 
density variation on the in-plane field, i.e., by coupling 
them through the Poisson equation) . The in- plane charge 
current Eq. H29|) and the the current leakage Eq. 125|) 
are to be solved jointly to specify completely the spin- 
dependent density profile N±{x) along the channel. We 
do this separately for the two systems under considera- 
tion. 



A. Silicon Inversion Layer 

A gate bias is necessary to create the inversion layer in 
silicon MOSFET-style device. We assume that the nor- 
mal metal gate is replaced with a ferromagnet, so that 
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the ferromagnet is biased with respect to the siHcon sub- 
strate. The Fermi level in the ferromagnet is lower than 
the Fermi level in the semiconductor. Hence 2DEG elec- 
trons can tunnel into the ferromagnetic gate, inelastically 
fall to the Fermi level in the ferromagnet, and have no 
way to get back into the semiconductor. This causes a 
leakage of the 2DEG density, so that the current is not 
constant along the channel, 



ox 



-3z,±{x, z = 0) = 



(31) 



Combining Eq. H31|l with Eq. H29(l . we have the differen- 
tial equation specifying the spin-dependent channel den- 
sity: 



d_ 

dx 



D±{x) 



dx 



eE,f±dN± N±ix) 



'Si,t 



dx 



. (32) 



The spin-dependent scattering times t± depend on the 
thickness of the oxide barrier under consideration. We 
have assumed that, except for at the source and drain, 
the system is homogeneous in x and that the lifetimes and 
scattering times are independent of the density variation 
throughout the channel. 

Because the diffusion constant is dependent upon the 
density in a complicated way, Eq. H32|l must be solved 
numerically. In the limit where the diffusion constant 
does not depend on x, the solution under a single gate is 

N±{x) = A±e-"/'±™ + B±e"/'± , (33) 
where the downstream and upstream length a^^'^'^ 



are 



down 
± 



2™si,t 

-eE^T±f± 
2™si,t 




{eE,YT±fl 
{eE,fT±fl 



+ 1 (34) 



1 (35) 



and A± and B± are fixed by the boundary conditions. 
Note that in all calculations in this paper, the in-plane 
electric field E^ was taken as negative, so that and 
^down (jgfj^ned to always be positive numbers. The 
real solution, which is calculated numerically with D±{x) 
varying along the channel, is more complicated, but the 
above approximation captures the qualitative aspects of 
the channel density. The density starts at its equilib- 
rium value near the source and the drain, which act as 
reservoir's for electrons. The electrons leak into the gate, 
causing the density to decrease as they move further from 
the source or drain. In the absence of a source-drain bias, 
the decay of the density would be symmetric about the 
center of the gate; a finite source-drain bias breaks the 
symmetry, and hence gives the two decay lengths and 
^down^ The spin-dependent channel density N±{x) deter- 
mines the spin-dependent channel current via Eq. (|29|l . 




1000 



400 600 
x-direction (Angstroms) 

FIG. 7: The spin-dependent 2DEG density (cm~^) and spin- 
dependent 2DEG current (cm~^ s~^) under a single ferromag- 
netic gate for the silicon inversion layer for zero source-drain 
bias (in-plane field of = V/cm). There is a diffusion 
current flowing from both the source and the drain into the 
gate. The other parameters used are explained in the text. 



The 2DEG density and current for the silicon inversion 
layer coupled to a ferromagnetic gate are shown in Fig. [7| 
and Fig. |Slfor the low- field regime E^ = V/cm and the 
high- field regime E^ — —5000 V/cm respectively. The 
source is at Xh = 0, the drain is at xr, = 1000 A, and 
both spin channels take their equilibrium density at the 
source and drain, N±{x'l) = N±{xji) = 0.5 ■ 10"'^^ cm~^. 
The oxide is assumed to be 6 A thick, so that the spin- 
dependent scattering times are t+ = 4 ps and r_ = 11 ps. 
The intrinsic scattering time is taken as 1 ps, and the 
temperature is 10 K. The spin-dependent density and 
current were calculated by solving the linear system of 
equations that results from the finite-differencing of the 
differential equation specifying the density fEa. l32ll . The 
solution was iterated until both the density N±{x) and 
the diffusion constant D±{x) converged. This method 
was used to calculate the current and density for the 
remainder of the paper. 

The effects of the leakage of channel carriers is im- 
mediately evident in Fig. [7| Carriers under the gate 
can irreversibly leak into the gate. A diffusion current 
flows from the source and from the drain into the 2DEG 
to replace the leaking electrons. The channel density is 
symmetric about the center of the device, and the cur- 
rent flows in opposite directions on the two sides of the 
gate. The spin-dependence due to the different scattering 
times for the two spin channels is evident; the -I- chan- 
nel, which leaks at a faster rate, falls to a lower density in 
the center of the gate, and has higher currents near the 
source and drain due to the larger gradient in the density 
as compared to the — channel. 

As the field is turned up, the drift term in Eq. H29|) 
fights against the backflow at the drain contact. At 
« —200 V/cm, the drift current overcomes the diffu- 
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The equation specifying the channel density is, from 
Eq. and Eq. 



200 400 600 800 

x-direction (Angstroms) 

FIG. 8: The spin-dependent 2DEG density (cm~^) and spin- 
dependent 2DEG current (cm~^ s~^) under a single ferromag- 
netic gate for the silicon inversion layer for an extremely high 
in-plane field of — —5000 V/cm. The other parameters 
used are explained in the text. 



sive backflow and the current through the drain becomes 
positive. Curiously, the field at which the drift current 
cancels the diffusion current is spin dependent, so that 
it is possible to create the following situations: (1) the 
drift current cancels the diffusion current for one of the 
spin channels, but not the other, so that the net current 
that flows through the drain is 100 % spin polarized, and 
(2) the current for one spin channel is positive while the 
current in the other spin channel is exactly opposite, so 
that no net current flows out of the drain, but a pure 
spin current flows out the drain. Further study of the 
pure spin current in the single ferromagnetic gate silicon 
system will be given in a future publication. 

In Fig. |S| a very high in-plane field is strong enough 
to overcome the backflow from the drain; the asymme- 
try induced by the strong source-drain bias is evident in 
the plot of the spin-dependent density. The transport is 
dominated by the drift term in Eq. H29|) due to the high 
in-plane held; the diffusion current is only a small correc- 
tion. More current is carried in the — spin channel than 
the -|- spin channel because the lifetime for the — spins, 
f_, is longer than the lifetime for the + spins, f-|_. 



B. InAs Surface Layer 

The 2DEG at the surface of InAs is natural and hence 
no gate bias on the ferromagnet is necessary. This keeps 
the Fermi levels in the semiconductor and in the ferro- 
magnet equal; no electrons will leak into the gate. The 
in-plane current in each spin channel must be conserved. 



dx \ dx J m' 



InAs 



dx 



(37) 



The spin-dependent density in the limit of D± (x) = con- 
stant is 



N±{x) = A±+B±e 



X I d± 



(38) 



dx 



'x,± 



(36) 



with A± and B± specified by the boundary conditions 
and the decay length d± = m\^p^^D±/ {—eExf±). The de- 
cay lengths are similar to the expressions found in Ref.fl^ 
and Ref. ^^for infinite spin lifetime. Because of the spa- 
tial dependence of the diffusion constant, the real solu- 
tion is more complicated, but the qualitative results still 
hold. 

The in-plane field drives a spin-accumulation at one 
of the boundaries. A boundary condition which speci- 
fies that the density is the same at the source and drain 
would prevent this from happening, resulting in a con- 
stant density throughout the channel. The diffusion term 
in Eq. H29(l vanishes, so that the current is only given by 
the drift term. The current carried by each spin chan- 
nel must be constant due to the absence of any 2DEG 
leakage. For zero source-drain bias, the current is thus 
zero in both spin channels. This is in marked contrast 
with the results for the silicon case. Fig. [7| in which a 
diffusion current flowed towards the center of the gate 
from both the source and the drain at zero source-drain 
bias. This simple behavior of the InAs surface layer with 
a single ferromagnetic gate is radically altered when two 
gates replace the single gate, as explained below. 



SPIN VALVE WITH TWO ADJACENT 
FERROMAGNETIC GATES 



Because of the general spin dependence in the 2DEG 
density and current as shown above, we now discuss a 
simple spin-valve type device to test the spin effects pre- 
dicted by the theory. The spin effects could not be seen 
in transport experiments with a single gate and nonmag- 
netic source and drain contacts. This is because the two 
spin channels are measured in parallel; switching the gate 
magnetization would exchange the roles of the two spin 
channels but would have no effect on the total current 
measured. The single-gate spin effects could in principle 
be measured through other spin-dependent means, but 
we concentrate here on transport effects that would be 
more useful in device considerations. 

To see the spin effects in transport experiments, we 
proposed a simple spin-valve design^^ shown in Fig. El 
The single gate in a normal MOSFET-style device is re- 
placed by two ferromagnetic gates that are very near each 
other. If a gate bias is applied (such as in the silicon sys- 
tem), the space between the two ferromagnets can be 
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FIG. 9: Schematic diagrams of the spin valve proposal in (a) 
the silicon and (b) the InAs systems. The device consists 
of normal nonmagnetic source and drain contacts and two 
ferromagnetic gates. A source drain bias Vsd creates the in- 
plane field Ex- In (a), a gate voltage is necessary to induce 
the inversion in the silicon system. The space between the 
ferromagnets is filled with a nonmagnetic metal. In (b), no 
gate voltage is necessary in the InAs system. An ultrathin 
oxide separates the ferromagnetic gates from the 2DEG in 
both systems, although the oxide can be omitted in (b). 



filled with a nonmagnetic metal such as aluminum in or- 
der to ensure that the 2DEG confinement is uniform be- 
tween the source and the drain. The nonmagnetic metal 
between the two ferromagnets is not necessary in sys- 
tems with no gate bias (such as the InAs system) . There 
is some anisotropy in the design of the two ferromagnets 
(either through geometry or material) so that the two fer- 
romagnets have different coercive fields, which allows the 
switching of the magnetization of the second ferromagnet 
while leaving the magnetization of the first unchanged. 
Because the two spin channels are effectively decoupled 
throughout the device (due to the long spin flip time), 
the addition of the second ferromagnet has a profound 
influence on the total current that flows in the 2DEG. 
This is because both the density and the current in each 
spin channel must be continuous throughout the device. 
This matching depends on the relative orientation of the 
magnetizations. The net effect is that the total current 
measured through the device will depend on the relative 
orientation of the magnetization of the two ferromagnets, 
which is a magnetoresistance effect. The first ferromag- 
net will be referred to as the first "gate" , and the second 
ferromagnet will be referred to as the second "gate" . 

Although we propose such a device as a way to test the 



predictions of the theory, the device as is could do per- 
form the role of a memory element. The magnetization 
of the first gate remains fixed, while the orientation of 
the second gate's magnetization depends on whether the 
information you are trying to read out is a "0" or a "1" . 
The measurement of the second gate's magnetization is 
read out by measuring the current at the drain contact. 
The bit could be written using a local magnetic field (cre- 
ated by nearby wires, as in current MRAM technology) 
that is strong enough to switch the second gate without 
switching the first gate. The two main benefits of such 
a design are that (1) the nonvolatile information storage 
has been incorporated onto the semiconductor, where the 
information processing occurs, and (2) we take advantage 
of the necessity for ultrathin oxide barriers due to the ag- 
gressive scaling of MOSFET technology to the nanometer 
scale (whereas currently the ability of 2DEG carriers to 
interact with the gate through an ultrathin oxide barrier 
is seen as a major obstacle to be avoided). Once demon- 
strated, other device designs can be explored which do 
not rely on analogy to existing structures to fully benefit 
from the new spin effects. 

Below we discuss the spin-valve for the silicon and 
InAs systems separately. To model the 2DEG density 
and current throughout the spin-valve, we assume the 
following. The source is at xl = and the drain is at 
— 2200 A. Both gates are 1000 A wide and the gap 
between the two gates is 200 A, so that the left side of 
the gap is at xa = 1000 A and the right side of the gap 
is at XB = 1200 A. 

The spin-dependent density with parallel gate magne- 
tization is N'^{x), where ± refers to the electron's spin 
under with respect to the first gate and p refers to paral- 
lel. The density and current must be continuous through- 
out the device. The spin-dependent scattering times are 
no longer constant along the channel due to the addition 
of the second gate, so that for parallel gate magnetiza- 
tions, 

{T± , XL < X < XA 
cx) , XA < 2^ < a;B , (39) 
T± , XB < X < X^ 

and the spin-dependent lifetimes are 

{f ± , Xl^ < X < XA 
To , XA < X <XB ■ (40) 
f ± , XB < X < Xji 

f± is specified by Eq. H28() . The density and current in 
both the silicon and InAs systems can now be calculated 
for parallel magnetization. 

For antiparallel gate magnetizations, the second gate 
magnetization is fiipped with respect to that of the first 
gate, so the spin-dependent scattering times for antipar- 
allel gate configuration are 

{T± , Xl, < X < XA 
oo , XA < X < xb , (41) 
Tzp , XB < X < xr 
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and the spin-dependent lifetimes are 



f± , Xh < X < XA 
To , XA < X < Xb 

f , Xb < X < xn 



(42) 



The obvious difference from the paraUel case is the ex- 
change of the roles of spin + and — under the second 
gate ( T± — » Tip and f± ^ fip). 

For both parallel and antiparallel gate magnetizations 
the density for both spin channels must take on its equi- 
librium value at the source and drain, which we assume 
for both the silicon and InAs systems is .5 • 10"'^^ cm~^ 
for each spin channel, so that the boundary conditions 
are 



Nlix^) 
Nlixn) 



(^l) 



0.5 • 10^2 cm-2 



(43) 



These boundary conditions are sufficient to fully calcu- 
late the spin-dependent density and current throughout 
the spin- valve for parallel and antiparallel gate magneti- 
zations. These are used next to calculate the density and 
current in the two-gate spin valve proposal for the silicon 
and InAs systems separately. 



A. Silicon Inversion Layer 

A schematic diagram of the silicon spin valve is shown 
in Fig. I^a). The addition of the second gate causes the 
scattering times and lifetimes to become dependent upon 
X. For parallel gate magnetizations, the differential equa- 
tion that must be solved is 



dx \ ^ dx 



. (44) 



"Si.t 



where the spin-dependent scattering time and lifetime 
along the channel are specified by Eq. H39() and Eq. H40() 
respectively. These equations must be solved numerically 
as previously described, with the boundary conditions 
Eq. gSJ. 

The density and current in the silicon inversion layer 
with parallel gate magnetization are shown in Fig. 1101 
The in-plane field that drives the current is — —500 
V/cm and the temperature is 10 K. The intrinsic scatter- 
ing time is To = 1 ps and the spin-dependent scattering 
times are t-^ — Aps and t_ = 11 ps, consistent with a 6 A 
barrier. The behavior is very different compared to the 
single gate case, Fig. |S1 The conductivity changes from 
under the first gate to the gap region (likewise from the 
gap region to under the second gate). Because E^ is con- 
stant, this requires that the density adjust itself to make 
both the density and current continuous. Spin accumula- 
tion can result if the in-plane driving field is sufficiently 
strong. Note that the current in the gap region of the 
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FIG. 10: The density (cm~^) and current (cm~^ s~^) for a 
silicon inversion layer with parallel gate magnetizations. The 
parameters are explained in the text. 
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FIG. 11: The density (cm~^) and current (cm~^ s~^) for a 
silicon inversion layer with antiparallel gate magnetizations. 
The two spin channels are labelled with respect to the first 
gate. 



device is constant because there is no electron leakage in 
this region. 

For antiparallel gate magnetizations, the differential 
equation that must be solved is the same as for the par- 
allel case (Eq. (Ejl ). except change p — > ap in all super- 
scripts. The scattering time is specified by Eq. (|41|) and 
the lifetime is specified by Eq. H42|) . Compared to the 
device with parallel magnetization, we see marked differ- 
ences in both the density and current in the antiparallel 
case (see Fig. Illf) . Again, this is because when crossing 
through the different regions of the device, the density 
and current must be continuous. Compared to the par- 
allel case, the -|- spin channel now sees the conductivity 
and leakage of the — spin channel under the second gate, 
so the matching is completely different. The total current 
is different in the two cases, causing a magnetoresistance 
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FIG. 12: The full line is the magnetoresistance (Eq. I|45^ ') as 
a function of the in-plane driving field E^c, using the same 
parameters as in Fig. 1101 and 1111 The dashed line is the 
magnetoresistance using scattering times appropriate for a 
thicker barrier. The high-field limit of this plot was discussed 
in Ref. 13 



effect. 

To clearly see the magnetoresistance effect caused by 
the switching of the second gate's magnetization, in 
Fig. El we plot the following magnetoresistance percent- 
age: 



MR = 100 



(45) 



where Jp(xr) and J!^^{x^) are the sum of the current 
in the two spin channels with parallel magnetization and 
antiparallel magnetization respectively. The current is 
assumed to be measured at the drain (x = xr). All 
parameters are the same as previously discussed. The 
magnetoresistance is quite small for large in-plane driv- 
ing fields, where the drift current dominates the trans- 
port. The difference between the parallel and antipar- 
allel currents, J|'(xr) — J^P(a;R,), is small compared to 
J^{xr) (or J^P(a;R)) because of the strong drift. As the 
in-plane field is decreased, the diffusion current becomes 
more important, which opposes the drift current at the 
drain. The different density profiles in the parallel and 
antiparallel cases imply that the diffusion current a the 
drain contact is different in the two cases; the spin effects 
are more pronounced and the magnetoresistance grows. 
At some critical field (w —200 V/cm) the drain current 
for parallel magnetizations equals zero; the drift current 
flowing in the -|-x-direction is exactly cancelled by the 
diffusion current flowing in the — x-direction. Due to the 
way in which Eq. 1)45(1 is defined, this causes a diver- 
gence of the magnetoresistance. As the in-plane field is 
reduced further, the total current in the parallel config- 
uration becomes negative as the diffusive backflow over- 
takes the drift current at the drain. At low source-drain 
bias the drain current in both the parallel and antiparal- 
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FIG. 13: The density (in units of cm~^) for an InAs surface 
layer with parallel gate magnetizations. The parameters are 
explained in the text. 



lei configurations are negative as the backfiow completely 
dominates the current; a small magnetoresistance is still 
present. The spin-valve would ideally be operated just 
above the critical field, where the drain current is still 
positive but the effects of the diffusion current are im- 
portant. In addition, the dashed line in the Fig.^lis the 
same calculation using the scattering times for a thicker 
oxide barrier. The effect is still present but moves to a 
lower in-plane field. This is because there is less leakage, 
and hence the diffusive backflow at the drain is lower 
than for a thinner barrier. 



B. InAs Surface Layer 

A schematic diagram of the InAs system is shown in 
Fig. Elb). For parallel magnetizations, the differenctial 
equation that must be solved is 



d_ 

dx 



^± dx 



eE^ d 
dx 



{rlNl) ^ , (46) 



'InAs 



where the scattering times and lifetimes are specified by 
Eq. H39|l and Eq. 140|l and the boundary conditions are 
specified by Eq. . 

The density in the InAs surface layer with parallel gate 
magnetization are shown in Fig. 1131 The in-plane field 
that drives the current and the boundary conditions on 
the source and drain densities are the same as in the 
silicon case discussed previously. The intrinsic scattering 
time is taken as 0.1 ps, and the spin-dependent scattering 
times are taken as t^. = 0.3 ps t_ =1 ps, consistent with 
a barrier of less than 5 A. In contrast to the single gate 
case, the density is not constant throughout the device, 
because it must adjust itself at the interfaces between 
the three regions to keep the current constant. Both spin 
channels must decrease their density in the gap region to 
keep the current constant, but the density profile for the 
two spin channels is very different because the lifetimes 
are different in the gate regions. This process leads to 
a static spin polarization at the interfaces between dif- 
ferent regions of the device, or spin accumulation. This 
accumulation must decay back to the equilibrium value 
at the drain. The current in each spin channel is constant 
throughout the device (1.47 ■ 10^* cm~^ s~^ for the -I- 
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FIG. 14: The density (cm~^) for an InAs surface layer with 
antiparallel gate magnetization. The two spin channels are 
labelled with respect to the first gate. The parameters are 
explained in the text. 

channel, 1.74 • 10^* cm^^ s^^ for the — channel) due to 
the "floating" gates. 

For the thinnest of oxides and intimate contact, Fig. El 
implies that the spin splitting becomes important. This 
would make the equilibrium values for the two spin chan- 
nels at the source and drain spin-dependent. This asym- 
metry further enhances the spin effects due just to the 
scattering time; the qualitative picture remains the same. 

For antiparallel gate magnetizations, the differential 
equation is the same as for the parallel case (Eq. 14t)|l 'l. 
except change p ^ ap in all superscripts. The scattering 
times and lifetimes are given by Eq. H41|l and Eq. (|42|l . 
As has already been discussed, the density and current 
are quite different from the parallel case due to different 
matching conditions for each of the spin channels. The 
density is plotted in Fig. Because of the exchanging 
of the roles of the two spin channels under the second 
gate, the density for the -I- channel increases to well over 
its equilibrium value in the gap region, and the — chan- 
nel decreases to well below its equilibrium value. Again, 
this is another example of spin accumulation at the inter- 
faces between different regions of the device. The current 
in each spin channel is constant throughout the device 
(« 1.6 • 10""^^ cm~^ s~^ for both spin channels). 

The magnetoresistance, Eq. (I45II . is relatively constant 
at « 1 % for all reasonable in-plane fields. This is because 
the diffusion term is always much smaller than the drift 
term (the spin accumulation that occurs is always very 
small, much smaller than the leakage-induced changes in 
the density in the silicon system) , so that the magnetore- 
sistance is basically just specified by the difference in the 
spin-dependent lifetimes. There is never any backflow, 
as in the silicon case, so the divergent structure in the 
silicon magnetoresistance (see Fig. 112(1 is not seen in the 
InAs case. 



VI. CONCLUSIONS 

We presented above a comprehensive theoretical treat- 
ment of the spin-dependent electronic and transport 
properties of a two-dimensional electron gas, under the 
strong infiuence of the proximity of a ferromagnetic layer. 



By constructing an appropriate Green's function, we de- 
termined the complex self-energy of the quantum con- 
fined electrons in the semiconductor, which are coupled 
quantum mechanically with the spin-polarized Fermi sea 
in the ferromagnet. Using a Hamiltonian coupling the 
two regions derived in a tight-binding-like approach, we 
calculated the spin-dependent properties of two paradig- 
matic systems: (i) the gate-induced inversion layer in a 
ferromagnetic metal-oxide-silicon junction; (ii) the (spon- 
taneous) accumulation layer of InAs separated from a 
ferromagnet by a thin oxide barrier. The ferromag- 
netic proximity induces a spin-splitting of the quan- 
tum confined subbands in the semiconductor and a spin- 
dependent broadening, which make the in-plane trans- 
port spin-dependent. We studied extensively the depen- 
dence of the ferromagnetic proximity as a function of the 
thickness of the thin oxide layer separating the semicon- 
ductor from the metal. Our results show that the spin- 
dependent lifetime broadening is the main effect, whereas 
the spin-splitting becomes sizeable only for nearly inti- 
mate contact between the semiconductor and the ferro- 
magnet. 

Knowledge of the spin-dependent electronic properties 
of the two-dimensional electron gas led to a treatment 
of the in-plane transport of spin-dependent current. The 
leakage current into the gate creates a density gradient 
along the semiconducting layer. The resultant drift and 
diffusion terms of the source to drain current in the semi- 
conductor were treated above in a self-consistent man- 
ner with the leakage and the density variation. While 
the leakage current is usually considered a limitation 
for electronic field-effect transistors with very thin ox- 
ide layers, the leakage current plays a positive role in the 
spin-dependence of the transport. We applied our trans- 
port theory to our recently proposed spin valve with two 
neighboring ferromagnetic gates. The detailed results of 
the spin-dependent steady-state densities and currents 
for different configurations of the magnetization in the 
gates yielded an explicit understanding of dependence of 
the magnetoresistance as a function of the source-drain 
bias. Notably, for a critical bias, a pure spin current 
(i.e., with a zero net charge current) can be created at 
the drain contact. This effect is caused by the gate leak- 
age (e.g. for the silicon inversion layer) and is due to the 
compensation of net drift and diffusion currents. 

We hope that the interesting new physics in the trans- 
port governed by the proximity effect will stimulate fur- 
ther explorations by experiments and by more realistic 
simulations, opening up the possibility of creating field- 
effect spintronics devices, as an alternative to spin injec- 
tion devices. 
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